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A magnetic field dependent Feshbach resonance has been 
observed in the elastic scattering collision rate between atoms 
in the F = 2, M = — 2 state of 85 Rb. Changing the magnetic 
field by several Gauss caused the collision rate to vary by a 
factor of 10 4 , and the sign of the scattering length could be 
reversed. The resonance peak is at 155.2(4) G and its width 
is 11.6(5) G. From these results we extract much improved 
values for the three quantities that characterize the interac- 
tion potential: the van der Waals coefficient C6, the singlet 
scattering length as, and the triplet scattering length ar- 



Very low-temperature collision phenomena can be 
quite sensitive to applied electromagnetic fields. Sev- 
eral groups have altered inelastic collision rates in op- 
tical traps by applying laser fields pL There have also 
been numerous proposals for using laser and static 
electric or magnetic fields to influence the S-wave scat- 
tering length (a), and equivalently, the elastic collision 
cross sections (a = 87ra 2 ) between cold atoms. Particu- 
larly notable is the prediction by Verhaar and co-workers 
H that as a function of magnetic field there should be 
Feshbach resonances in collisions between cold (~ [iK) 
alkali atoms. These resonances were predicted to cause 
dramatic changes in the cross section, and to even al- 
low the sign of the scattering length to be changed. Such 
resonances are very interesting collision physics, but they 
also offer a means to manipulate Bose-Einstein conden- 
sates (BEC), the properties of which are primarily de- 
termined by the scattering length. Some time ago we 
searched for such resonances in 133 Cs and 87 Rb, without 
success This was not surprising because there was 

enormous uncertainty as to the positions and widths of 
the predicted resonances. However, recent photoassoci- 
ation spectroscopy has greatly improved the knowledge 
of the alkali interatomic interaction potentials, allowing 
these resonance parameters to be predicted with much 
less uncertainty . As we show below, the fact that the 
widths and positions of these resonances are so sensitive 
to the interatomic potentials makes their measurement a 
good way to determine these potentials. 

The improved predictions for the resonance positions 
facilitated their observation. In the past few months Fes- 
hbach resonances have been observed in both 23 Na ]Tc| ] 
and 85 Rb p[ . In the sodium work, two magnetic field 
dependent processes were observed: 1) a change in the 



expansion rate of BEC due to a change in the scattering 
length, and 2) an enormous, and as yet unexplained, in- 
crease in the loss rate. This loss rate precluded the study 
of collisions in the interesting field regime near the cen- 
ter of the resonance where the scattering length changes 
sign; the width of the resonance was not measured. In 
the rubidium work [Til ] the resonance was detected as a 
magnetic field dependence of the photoassociation spec- 
tra. The resulting resonance width was measured to be 
far larger than originally predicted. 

Here we report the study of a Feshbach resonance in 
the elastic collision cross section between atoms in the 
F = 2, M = -2 state of 85 Rb. By changing the DC 
magnetic field we are able to change the collision rate by 
4 orders of magnitude, and explore regions of positive, 
negative, and essentially zero scattering length. We de- 
termine the width and position of the resonance about a 
factor of ten more precisely than in Ref. and from 
this data we improve the accuracy of the Rb interaction 
potential substantially. In contrast to Ref. pfj[ , we do 
not observe two or three body loss at resonance because 
we work at much lower densities (10 9 atoms/cm 3 ). 

We measured the collision rates using the technique of 
"cross dimensional mixing," as in our earlier work [Q. 
In this technique, a nonisotropic distribution of energy 
is created in a magnetically trapped cloud of atoms, and 
the time for the cloud to reequilibrate by elastic collisions 
is measured. The apparatus used is identical to that used 
in our previous work on BEC in 87 Rb p^ |. It is a dou- 
ble magneto-optic trap (MOT) system in which multiple 
samples of atoms are trapped in a relatively high-pressure 
chamber, and then transferred to a second MOT in a low- 
pressure chamber. The second MOT is then turned off 
and a "baseball coil" magnetic trap is turned on around 
them. The atoms are then cooled by forced evaporation. 

After the atoms are evaporatively cooled to the de- 
sired temperature, the energy in the radial direction (and 
correspondingly the square of the width of the trapped 
cloud) is reduced to about 0.6 that of the axial direc- 
tion. This is done by decreasing the frequency of the 
RF "knife" used for the evaporative cooling more rapidly 
than the cloud can equilibrate. The bias magnetic field 
is then adiabatically ramped to the selected value, and 
the cloud is allowed to equilibrate for a fixed time. The 
shape of the cloud is then measured by absorptive imag- 
ing. This is repeated for different equilibration times, 
and the aspect ratios of the cloud vs. time are fit to 
an exponential to determine the equilibration time con- 
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stant r eq . As discussed in Ref. O], r eq is proportional 
to the inverse of the elastic collision rate . From the 
measured trap oscillation frequencies and the measured 
shape and optical depth, we determine the temperature T 
and average density (n) of the cloud. The magnetic field 
value at the center of the cloud is determined to ±0.1 G 
by finding the RF frequency that resonantly drives spin 
flip transitions for the atoms at that position and using 
the Breit-Rabi equation. The spread in magnetic field 
across the cloud scales as T 1/2 , and is 0.6 G FWHM for 
a 0.5 fiK cloud. 

The measured equilibration times depend on B and 
T and vary from 0.15 s to nearly 2000 s. To display 
data taken at several values of both density and temper- 
ature, we convert it to the normalized equilibration rate, 
r norm = l/{nv)T eq . Figure 1 displays the data versus 
mean magnetic field. T norm has units of cross-section 
and in fact is proportional to an empirical average over 
the field- and energy-dependent elastic scattering cross- 
section a(B,E). To calculate this average is problem- 
atic — a detailed Monte Carlo simulation of particle tra- 
jectories JlHj is probably required. However in the low- 
temperature limit it is clear that the maximum and min- 
imum values of T norm will occur at the same magnetic 
field as the peak and zero of a(B,E = 0). For a rough 
qualitative comparison, in Fig. 1 we have also plotted the 
theoretically predicted values of a(B, E) for two temper- 
atures. These predictions are calculated using the atomic 
potentials determined from the observed maximum and 
minimum values of the relaxation rate, as discussed be- 
low. Both the cross section curves and the equilibration 
data show the same qualitative features. For temper- 
atures below a few /xK, there is a slightly asymmetric 
peak near 155 G, and the width and height of this peak 
are strong functions of temperature. At 167 G there is 
a profound drop in the rate. This dip is quite asymmet- 
ric, but the shape is relatively insensitive to T, and at 
the bottom (field value B m i n ) the rate is essentially zero. 
The field value of the peak (-B pca k) is customarily defined 
to be the position of the Feshbach resonance, and the res- 
onance width, A is the distance between -B poa k and B m i n 
flTf . The scattering length is positive for field values be- 
tween Speak and Bmin, and is negative for field values 
below -Bpcak or above -B m ; n . This dependence of the sign 
is expected from previous theory. The observation that 
-Speak is at a lower field than -B m ; n provides experimental 
confirmation that the scattering length is negative away 
from the resonance. 

Because the functional form of the normalized equili- 
bration rate is not known, the desired quantities, -B pca k 
and Bmin cannot be found by a detailed fit to all of the 
data. We determine them by fitting a simple smooth 
curve to only the few highest (or lowest) points at each 
temperature below 5 /xK, and assigning correspondingly 
conservative error bars that more than span the values 
determined for all temperatures |Q. We find B pGa k = 



155.2 (4) G and -B m ; n = 166.8 (3) G, giving a resonance 
width A = 11.6 (5) G. The values for -B pea k and A are 
reasonably consistent with the less precise values 164(7) 
and 8.4 (3.7) G measured in Ref. pj. In our experi- 
ment, the accuracy of the peak position is much higher 
primarily because of better field calibration. The better 
accuracy in the width is largely because the photoasso- 
ciation measurement is a somewhat less direct way to 
observe the resonance, and so substantial and somewhat 
uncertain corrections are required to go from observed 
signal widths to the actual resonance width. 

We can now use these measured quantities to deter- 
mine the singlet and triplet Born-Oppenheimer poten- 
tials. The accuracy of predicted cold collision observables 
hinges on the quality of these potentials used in the radial 
Schrodinger equation for the nuclear motion. Accurate 
determination of these potentials in turn relies primarily 
on three parameters, the long range van der Waals coef- 
ficient Cg and the zero-energy singlet and triplet phases. 
The phases are usually tabulated in terms of singlet as 
and triplet ax scattering lengths calculated with hyper- 
fine terms omitted from the radial equation. 

The Feshbach resonance involves states with both sin- 
glet and triplet character; consequently the present de- 
termination of resonance width and position versus B 
can be translated with the aid of theory into greatly 
improved values for the singlet and triplet scattering 
lengths. The old "nominal" values of the scattering 
lengths and Cg (optimized to achieve agreement with 
previous measurements) predicted the position of this 
Feshbach resonance reasonably well, but, as in Ref. p"l| ] 
the predicted width was much smaller than observed. 
This width reflects the coupling of bound and contin- 
uum channels, and is primarily controlled by the differ- 
ence between the singlet and triplet scattering lengths. 
It is inconvenient to tabulate results as a function of 
singlet and triplet scattering lengths for 85 Rb, because 
both are very close to lying on top of a divergent pole. 
Accordingly, we present the results of our analysis in 
terms of a dissociation phase vd as in Ref. @. We 
define the corresponding singlet i^ds or triplet i^dt dis- 
sociation phase in terms of the scattering length through 
the relation [j^: i^d = (1/ 7r) cot _1 (a/a re f — 1) where 
a rof = r(5/4)/[V8r(3/4)](2mC 6 /n 2 ) 1 / 4 g| and m is 
the reduced mass of the atomic pair. The dissociation 
phase is related to the short range quantum defect /i sr 
presented in Ref. |ll| by z/d = /i sr — 1/4. The quantum 
defect method developed in Ref. p9] was used to generate 
the theory curves in Fig. 1. 

We have adjusted the singlet and triplet Rb-Rb poten- 
tial curves at small distances until the theoretical calcu- 
lation agrees simultaneously with the present measure- 
ments of -Bp Ca k and A, to within the experimental un- 
certainties. For a given Cq, this severely constrains the 
values of both z^ds and z^dt (see Fig. 2). However, the 
position and width are not sufficient to totally constrain 
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all three parameters C6, vds, and vdt- In order to ob- 
tain values for all the parameters, we also require that 
our triplet potentials predict the g-wave shape resonance 
in an energy range consistent with the measured value 
given in Rcf. |2C[ . 

The results of this analysis are presented in Fig. 2. It 
is immediately clear from the figure that accurate mea- 
surement of Feshbach resonances are an extremely precise 
method for determining threshold properties of the inter- 
atomic potentials. In fact, our measurement has allowed 
us to reduce the combined Gq, vns, and z^dt parameter 
space by roughly a factor of 80. As discussed in Ref. [ Q |, 
the position of the resonance depends mostly on C6 and 
the sum ^ds + ^dt- The width depends mostly on the 
difference vbs ~ ^dt- As Fig. 2 shows, the allowed vus, 
i»dt parameter region is an extremely correlated function 
of Cq. In particular, we find that the nominal value of 
the sum and difference coordinates are linear functions of 
Cq and the uncertainties in the coordinates are constant. 
Accordingly, we can represent the allowed parameter re- 
gion in the following manner: 

J^ds + ^dt = -0.0442 + 2.2(10~ 4 )(C 6 - C 6 ) ± 0.0005 and 
vps - fDT = 0.0666 - 6.2(10- 5 )(C 6 - C 6 ) ± 0.003 where 
C6 = 4700. Converting the dissociation phases into scat- 
tering lengths, we find reasonable agreement with previ- 
ous work [^1-23|. We find for the nominal C6 the val- 
ues (in a.u.) a t ( 85 Rb)=-363±10, a s ( 85 Rb)=2300±^, 
a t ( 87 Rb)=109.3±0.4, and a s ( 87 Rb)=92.7±0.4. Our value 
for Co, 4700(50) a.u. is slightly higher (and with smaller 
uncertainty) than a previous analysis based solely on 
the <?-wave shape resonance |2^|. We have confirmed 
that several scattering observables predicted by the new 
Born-Oppenheimer potentials are consistent with previ- 
ous measurements. Specifically, the new potentials pre- 
dict a broad d-wave shape resonance 
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scattering length ratio ct2,i/fli,-i |26J] and the thermally- 
averaged |2, 2) + |l,— 1) inelastic rate constant ^71 that 
are consistent with previous measurements. We also find 
10 of the 12 measured d-wave bound states |22| within 
the 2a error bars The new potentials also permit 

us to predict additional 85 Rb Feshbach resonances ||, at 
S pcak = 226.5 ± 4.8 G with a width of 0.01G, and at 
S pcak = 535.8 ± 4.0 G with a width of 2.2 ±0.2 G. 

We have shown that precise measurements of the Fesh- 
bach resonance in the elastic scattering cross section pro- 
vide unprecedented knowledge of the Born-Oppenheimer 
potentials that control Rb-Rb scattering processes at 
sub-mK temperatures. We have also used this reso- 
nance to enhance evaporative cooling of 85 Rb in a low 
density magnetic trap, and change the scattering length 
from negative to positive. These capabilities should allow 
novel studies of BEC in the future. 
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FIG. 1. (a) The data points show the measured equi- 
libration rate divided by average density and velocity 
(F„. orm = l/{nv)r eq ) vs. magnetic field B. Data are shown 
for clouds at five different temperatures: 0.3 /iK (filled tri- 
angles), 0.5 /iK (•), 1.0 /iK (o), 3.5 /iK (v), and 9.0 /iK 
(filled squares). The lines are calculations for the thermally 
averaged elastic cross sections (not equilibration rates!) and 
hence are not expected to fit the data points. The solid line 
corresponds to 0.5 /iK and the dot-dashed to 9.0 /iK. (b),(c) 
Expanded views of the regions of maximum and minimum 
cross section. 

FIG. 2. Comparison of the allowed uds, vdt parameter 
space based on this work and previous measurements for a 
constant value of C6- The rectangle is the allowed range from 
Ref. [22]. The large diamond is the range from Ref. [11]. 
The small filled diamond is the range from this work. Note 
that the rectangle and large diamond constraints were estab- 
lished using a different value of C6 from the one used in this 
work (4550 instead of 4700 a.u.). As is discussed in the text, 
the position of the constrained regions depends linearly upon 
C6, which is indicated by the small connected diamonds on 
either side of the small filled diamond. These diamonds show 
the effect of 50 a.u. uncertainty in C6- The arrows indicate 
sensitivity to changes in the resonance width and position in 
this paramter space. 
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